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We propose a hybrid quantum repeater protocol combining the advantages of continuous and 
discrete variables. The repeater is based on the previous work of Brask et al. [Phys. Rev. Lett. 
105, 160501 (2010)] but we present two ways of improving this protocol. In the previous protocol 
entangled single-photon states are produced and grown into superpositions of coherent states, known 
as two-mode cat states. The entanglement is then distributed using homodyne detection. To 
improve the protocol, we replace the time-consuming non-local growth of cat states with local 
growth of single-mode cat states, eliminating the need for classical communication during growth. 
Entanglement is generated in subsequent connection processes. Furthermore the growth procedure is 
optimized. We review the main elements of the original protocol and present the two modifications. 
Finally the two protocols are compared and the modified protocol is shown to perform significantly 
better than the original protocol. 

PACS numbers: 03.67.Bg, 03.67.Hk, 42.50.Ex 



I. INTRODUCTION 

A major goal in the field of quantum information is 
distributing entanglement over large distances. A strong 
motivation for this is that it may enable transmission of 
information, which is secure against eavesdropping [TJ [2] 
even in cases where the measurement devices or the 
source are untrusted [3]. More generally the distribu- 
tion of entanglement is required for almost any task in 
quantum communication. Direct distribution of entan- 
glement requires transmission of fragile quantum states, 
which is difficult in practice due to loss and decoher- 
ence in optical fibers. Quantum repeaters overcome this 
problem by initially generating entanglement over short 
distances and then distributing it via entanglement swap- 
ping, which only requires local operations and classical 
communication [4]. Recently, much effort has been de- 
voted to the construction of quantum repeaters based 
on atomic ensembles [5HI9], Despite the experimental 
advances towards this goal the construction of a fully 
functioning quantum repeater remain a formidable chal- 
lenge. In particular the low efficiencies obtained in cur- 
rent experiments severely reduce the achievable commu- 
nication rates. Numerous theoretical proposals have been 
developed to improve the communication rates [7][T8ll20l- 
[22] . In particular people have suggested schemes relying 
on multiplexing of operations so that high local repeti- 
tion rates can compensate for the slow non-local oper- 
ation which require communication with distant parties 
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Quantum communication generally works in two 
regimes; the discrete and the continuous variable regime. 
In the discrete variable regime, information is carried by 
single photons and measurements rely on single photon 
detection (SPD). This facilitates detection and correction 
for loss, but the efficiency of most available single-photon 
detectors is low, reducing the rate of entanglement distri- 
bution. High-efficiency (> 90%) SPD is possible, but re- 



quires detectors, such as superconducting transition-edge 
sensors, which are expensive and not widely available 
[25] . In the continuous variable regime, information is 
encoded in operators with a continuous spectrum such as 
the field quadratures of the electromagnetic field. These 
are measured using homodyne detection which is very ef- 
ficient in practice (~ 99%) but has the drawback that loss 
is not as easily detected as in the discrete variable regime. 
Recently hybrid quantum repeater protocols, combining 
the two regimes, were proposed first for spin systems in 
cavities [26] and later for atomic ensembles [27]. Here 
we follow the approach of Ref. [27 for atomic ensembles. 
The performance of this repeater protocol is comparable 
to the best proposed atomic-ensemble based repeaters in 
the discrete variable regime if these are operated using 
realistic SPD with limited efficiency [5j [17] [18] . 

Here we propose and analyze two modifications to the 
protocol of Ref. [27] . The hybrid repeater protocol cre- 
ates entanglement between two stations in the form of 
single-photon superpositions 1 01) + 1 10) ignoring normal- 
ization for simplicity. Through a probabilistic procedure, 
these states are then grown into states resembling 

| 7 (0, a)) cx e iG \a) a \a) h + e~ iG \ - a) a \ - a) b , (1) 

by means of local operations and classical communica- 
tion. Here \a) a denotes a coherent state with amplitude a 
in mode a, and is a phase. We refer to states of the form 
in Eq. as two-mode cat states since they are two mode 
superpositions of two "classical" states |a), | — a). Be- 
cause classical communication between distant stations 
is time consuming, the growth procedure is slow. In par- 
ticular the single-photon entanglement generation step, 
which has low success probability, needs to be repeated 
every time the growth step fails. In a related setup a 
solution to this problem was suggested in Ref. [15]. To 
improve the communication rate, it was proposed to re- 
place the low success non-local entanglement generation 
by a rapid preparation of a suitable local states. Because 



local operations do not rely on communication with dis- 
tant parties they have a much higher obtainable rate. 
The locally generated states are more suitable for entan- 
glement generation and can be connected with a much 
higher probability reducing the time spent on the slow 
non-local operations. Here we follow a similar path and 
consider interchanging the first two steps of the repeater 
protocol such that states resembling one-mode cat states, 

\£(e,a))cxe i0 \a)+e- i9 \-a), (2) 

are first grown locally and then subsequently connected 
to create entanglement by means of non-local single- 
photon subtraction. Such a modification reduces the 
need for classical communication and allows a higher rep- 
etition rate to be reached. This is the main idea behind 
the new repeater protocol detailed below. In addition, we 
optimize the cat-state growth procedure, improving the 
rate further. The remainder of the paper is organised 
as follows. In Sec. [IT] we review the scheme of Ref . [27] . 
In Sec. |III| we present the modified protocol. In Sec. |IV| 
we present results of numerical simulation of the new 
scheme, evaluating the performance in terms of achiev- 
able rates. Finally we conclude in Sec. (V) Some detailed 
calculations and parameters relevant to the numerics are 
given in appendices [Aj [Bj [C] and |DJ 

II. REVIEW OF PREVIOUS SCHEME 

The repeater protocol of Ref. [27 consists of three 
steps, (i) heralded entanglement generation based on 
sources of two- mode squeezed vacuum and SPD, (ii) 
growth of two-mode cat states from entangled single pho- 
tons by means of homodyning, and (iii) entanglement 
swapping based on homodyning. The steps are outlined 
in Fig. [I] 

In step (i) (see Fig. |TJi) ) , two sources produce two- 
mode squeezed vacuum states of the form 

|0,0) + V&IM) +0(p), (3) 

where p pa ir is the probability to produce a photon pair. 
These sources can be realized using parametric downcon- 
version crystals or ensembles of A- type atoms [5j[l6]. One 
output mode from each source is read into a quantum 
memory while the remaining modes are sent to a bal- 
anced beam splitter positioned between the two sources. 
The beam-splitter outputs are measured and a single 
SPD click projects the two modes in the quantum mem- 
ories into an entangled state |0, 1) + |1,0). The pair- 
production probability p pa ir (and hence the squeezing) 
needs to be small to ensure that the final state does not 
contain more than a single photon. 

In step (ii) (see Fig.[TJii)), two entangled single-photon 
states are combined on balanced beam splitters and the 
X quadratures of one output mode from each beam split- 
ter are measured (a related procedure to perform distil- 
lation of continuous-variable entanglement was shown in 




FIG. 1. (Colour online) Steps of the protocol in Ref. [27]. 
(i) Entanglement is generated using two sources of two-mode 
squeezed vacuum. One mode from each source is transmitted 
to a balanced beam splitter and the outputs are measured. 
Detection of a single photon heralds entanglement between 
the remaining modes stored in quantum memories (QM). (ii) 
Growth of cat states. Two entangled states are combined 
locally on balanced beam splitters and the X quadrature is 
measured. Success is conditioned on the sum of the outcomes 
taking a value close to zero, (iii) Entanglement swapping. 
One mode from each state is combined on a balanced beam 
splitter and the X and P quadrature of the outputs are mea- 
sured. Success is conditioned on a value of the A-outcome 
close to zero. 

Ref. [28j [29]). Whenever the sum of the measurement 
outcomes fulfill \x a + < A, for a certain acceptance 
interval A, the state is kept. The process can be iterated 
by combining two states resulting from successful growth 
at the previous level and repeating the procedure. In the 
limit of small A, the final output resulting from this pro- 
cedure approaches a non-locally squeezed two-mode cat 
state of the form 

5+(2)| 7 (0,M m /v / 2)>, (4) 

where \i m = ^2 m + 1/2 and m is the number of itera- 
tions. 5+ (2) denotes non-local squeezing in the variance 
of X a + Xb by a factor of two. The squeezing operator 
has the general form S(() = exp [\C^ 2 ~ \C<^ 2 )- The 
acceptance interval A determines the probability for suc- 
cessful growth, and hence the rate, as well as the fidelity 
of the output state with respect to the state in Eq. Q. 
Larger A corresponds to higher success probability but 
lower fidelity. The choice of A thus defines a tradeoff be- 
tween the rate and the fidelity. In Ref. [27 , A was fixed 
to take the same value for all m. 

The final step (iii) (see Fig. []Jiii)) is entanglement 
swapping where neighboring entangled segments are 
combined to create longer segments. Two modes, one 



from each entangled pair, are combined on a balanced 
beam splitter and the X and P quadratures of the out- 
put modes are measured. Whenever \x\ < 5, x being 
the outcome of the X measurement, the entanglement 
swapping is considered to be a success and the output 
state is kept. The process is iterated until entanglement 
is distributed over the total length, L of the repeater. 
This is obtained by first dividing L into 2 n segments of 
length Lo = L/2 n over which entanglement is created. 
At each swap level, every two neighbouring segments are 
connected, such that the entanglement distance is dou- 
bled. After n swap levels entanglement is distributed over 
the entire length L. In the limit of small J, the state pro- 
duced after n swap levels approaches a locally squeezed 
two-mode cat state 



\Tpideal)ab = S a (V2)S b (V2) | 7 (</>„, l^l* ^m)) ab, (5) 



where the phase (j) n depends on the P-measurement out- 
comes from the previous levels. This state contain one 
ebit of entanglement and is used to quantify the perfor- 
mance of the repeater via the fidelity 



(ii) 




FIG. 2. (Colour online) Steps of the modified repeater, (i) 
Local growth of cat states. The output modes from two 
sources of two mode squeezed vacuum states are combined 
on a balanced beamsplitter and the X quadrature of one of 
the outputs is measured conditioned on a click in both SPD- 
detectors. Conditioned on the measurement outcome the re- 
sulting state is kept for further processing, (ii) Entangle- 
ment generation from single-mode cat states. Small parts are 
tapped off from two input cat states at two separate loca- 
tions, and the remaining parts are stored in quantum memo- 
ries (black dots) . The fraction tapped off is controlled by the 
reflectivity r of the local beam splitters. The reflected signals 
are transmitted to a central, balanced beam splitter, and the 
output ports are measured. Conditioned on a click in either 
of the detectors, the memories are prepared in an entangled 
state. 



F pr ev = \(^ideal\p\^ideal)\ , (6) 

A. Growth of cat states 



where p is the density matrix of the final output state 
of the repeater. As for the growth step (ii), there is a 
tradeoff between fidelity and rate through the acceptance 
parameter S. The upper limit to the success probability 
of entanglement swapping is 1/2 for the simple proce- 
dure considered here but the success probability can be 
increased by using a more complicated procedure. For 
simplicity we will not consider this here. 



III. THE MODIFIED SCHEME 



To improve the rate of entanglement distribution, we 
interchange steps (i) and (ii) above, resulting in a new 
protocol based on local growth of single-mode cat states 
and subsequent non-local single-photon subtraction. The 
steps of the new protocol are sketched in Fig. [2] Entan- 
glement swapping is performed in the same manner as in 
step (iii) above. 



The first step of the modified protocol is growth of 
states approximating squeezed single-mode cat states 

\( m ) = S(2)— = (\^ m ) + \-^ m )), (7) 

where S(2) denotes squeezing by a factor of two in the 
variance of the X quadrature, |/i m ) is a coherent state 
with amplitude /i m = y/2 m + 1/2 and N+ m is a nor- 
malization constant. These states can be grown by a 
setup very similar to step (ii) of the original protocol, 
as explained in Ref. [27 . The input states are single- 
mode one-photon states, generated by detecting one half 
of a two-mode squeezed state with small pair-production 
probability p. To understand the growth procedure, we 
consider the ideal limit where each source produces a pure 
single-photon state |1) with corresponding wave function 
for the ^-quadrature 

^oW = ^^e'^x. (8) 

The joint wave function before the beam splitter is 
ipo(x)ipo(y). At the output of the balanced beam splitter, 
this is transformed into ip^x + y) / \/2)ij)o((x — y)/y/2) oc 
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e 2<> 2 +^ 2 )(£ 2 — y 2 y Now mode y is measured and the is approximated by 



state is kept if yo G [—A, A] where yo is the measure- 
ment outcome. Taking the limit A — >> we find the out- 
put state ^i(x) cc e~z x x 2 . Then the process is iterated 
with tpi(x) as input. After m iterations the output wave 
function becomes 



(9) 



The overlap of this state with the state in Eq. ^ exceeds 
99% for m > 2 and approximate squeezed cat states can 
thus be grown this way. 

As in the previous section, there will be a tradeoff 
between the fidelity and the rate controlled by A. In 
Ref. [27 A was kept fixed at the same value in every 
iteration but here we investigate the improvement by al- 
lowing different values of A for each m. To understand 
the possible improvement allowed by varying the interval 
we first analyze how the growth procedure works. The 
output wave function of the growth procedure (approxi- 
mately ijjm{x) for small A) is symmetric with two peaks; 
one at x < and one at x > 0. Suppose that the mea- 
surement is performed in the symmetric output of the 
beam splitter with quadrature operator X + = X\ + X2 • 
In this mode the quadratures add. If the two peaks with 
positive x are combined the measurement outcome will 
likely have a positive value. Similarly, combining the 
negative peaks leads to a negative outcome. These two 
possibilities are not desirable since the wave function in 
the antisymmetric mode X_ = X\ — X2 essentially will 
be a peak around zero, because in this mode the quadra- 
tures subtract. However, when a negative and a positive 
peak combine, the measurement outcome will be in the 
vicinity of zero. Since there are two paths leading to 
this result, corresponding to two different states in the 
antisymmetric output mode, the desired cat state is gen- 
erated. The acceptance interval must be chosen such 
that one avoids outcomes resulting from the tail of the 
distribution coming from the combination of two posi- 
tive or two negative peaks. The closer the peaks are to 
each other at the input, the smaller acceptance interval is 
allowed. As the growth process is iterated, the peaks be- 
come more separated and larger acceptance intervals can 
be chosen, resulting in a higher probability of success. 

We have optimized the acceptance interval to achieve 
the highest possible probability for a fixed target fidelity 



-Fgrowth — (Cm I Pm \Cm) 



(10) 



of the output state p m of the growth procedure. We 
assume perfect one-photon states at the inputs, and 
calculate the fidelity and rate on a grid of values for 
each acceptance interval, A m under the constraint that 
A m +i > A m . The optimization was made using Wigner 
functions, since these provide a natural description of 
mixed continuous- variable states and make it possible to 
compute the average output fidelity. Details are given in 
App. [Al The rate, in units of the source repetition rate 



growth 



m—1 



PiP 2 -..P m , 



(ii) 



where P m is the probability of successful growth in iter- 
ation m. This expression assumes that the outcomes of 
successful events can be stored while unsuccessful events 
are repeated until they succeed [18 . Note that by assum- 
ing that successful events are stored in quantum mem- 
ories we avoid the usual exponential scaling with the 
number of conversion events (recall that the number of 
down conversions is 2 m ). The result of the optimization 
is shown in Fig. [3] where we plot F growt h against Rgrowth- 
Note that the fidelity does not reach unity as R gr0 wth — > 
since we take the fidelity with the approximate cat state 
in Eq. ^ and not the state in the ideal limit in Eq. 
([9]). The calculation was restricted to m < 5 for runtime 
reasons. 
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FIG. 3. (Colour online) Optimized production rate of ap- 
proximate squeezed cat states after (a) m=2,3 iterations (b) 
m=4,5 iterations. The dashed lines are the optimal curves 
and the solid curves are obtained for identical acceptance in- 
tervals in all iterations. The fidelity was calculated with the 
target state in Eq. |7| and the rate is given in units of the 
source repetition rate. 
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Fig. [3] shows that the growth procedure is indeed im- 
proved by allowing for different acceptance intervals in 
every iteration. However the rate is not significantly im- 
proved for a small number of iterations m. Nonetheless 
Fig. [3] indicates that the improvement will increase with 
m since for larger m the peaks in the input states become 
more separated. For an output fidelity of 0.9 the ratio of 
the modified rate to the previous rate is 1.03, 1.10, 1.21, 
and 1.53 for m = 2,3,4, and 5 respectively. For very 
large m the two peaks will be so far separated that we 
can choose an acceptance interval for which the success 
probability approaches \ without affecting the fidelity. 



B. Connection of cat states 

The second step of the new protocol is to create entan- 
glement by connecting the single-mode cat states from 
the first step. The method we employ was proposed by 
N. Sangouard and coworkers in Ref. [17] and the setup is 
shown in Fig. [2|ii). 

A small part is subtracted from each input state by 
means of asymmetric beam splitters with low reflectiv- 
ity r. The remaining parts of the states are stored and 
the reflected parts are transmitted to a central station. 
Here, the two signals are combined on a balanced beam 
splitter and the two output ports are measured with 
photodetectors. Successful entanglement generation is 
conditioned on a click in exactly one of the detectors. 
When a click is observed in the symmetric output port 
of the beam splitter, the quantum memories are projected 
into an entangled state approximating the two-mode cat 
|7(0, ay/1 — r)). For a click at the antisymmetric output, 
the state is identical up to a local phase shift. The pro- 
cedure can be understood easily in the ideal case where 
the inputs are exact cat states |£(0, a)) and r — » 0. In 
this limit, a click heralds non-local subtraction of a sin- 
gle photon from the joint state of the memories. The 
memories are thereby projected into the (unnormalized) 
state 



(a±6)|e(0,a)) a |e(0,a)) 6 , 



(12) 



where a, b label the output modes. Inserting the defini- 
tion ([2| and recalling that coherent states are eigenstates 
of the annihilation operators, we notice that the compo- 
nent of the wave function where the a and b mode have 
the opposite (same) phase e.g. \a) a \ ~ a )b (\ a )a\ a )b) 
vanish by interference for the plus (minus) combination. 
Therefore the resulting state is 

\ot) a \ ±Ot)b-\- Ot) a \ =F Ot) b . (13) 

Comparing to Eq. ([I]), we see that up to a local phase 
shift this state is equal to \^(tt/2, a)). 

In practice, it is very hard to create genuine cat states. 
Therefore we shall use the approximate squeezed cat 
states from the previous step of the repeater protocol. 
We now examine the behaviour of these states under en- 
tanglement generation in the ideal limit A — >• 0, in which 



case they are given by Eq. ([9| , and taking again the limit 
r — » 0, we find 

\^m)ab = (CL± b)\^m)a\^m)b, 

1 (14) 

(X —r(\0m) a \lm)b ± |lm)a|0 m )&), 

V2 



where 



<x|0 ro ) = (x\ a = T(2 m - l/2)- 1 /V m - 1 e-^ 2 
<x|l ro ) = t/> m (x) = T(2 m + 1/2)" 1 / W 



(15) 



Here, |l m ) is a superposition of even photon states and 
|0 m ) is a superposition of odd photon states, and for m > 
2 they resemble squeezed even and odd single-mode cat 
states respectively i.e: 



|lm>«£(2) 



|0 m ) « 5(2) 



Mm)) 



(16) 



(17) 



where jl ri 



y/ 2m ~ 1A Vri 



y/2 m + 1/2 and 



AT+ , N~ are normalization constants. For m > 2 the 
fidelities between |l m ), |0 m ) and the respective cat states 
are both > 99%. 

The state |\P m ) contains one ebit of entanglement and 
is obtained in the low-rate limit of small acceptance inter- 
vals during growth and small reflectance during connec- 
tion. |\I/ m ) however deviates from the squeezed two- mode 
cat state that was shown to be useful for entanglement 
swapping in Ref. [27 . For m = 2 the overlap with a 
locally squeezed two-mode cat state of the form used in 
Ref. [27] 



£(2) a £(2) 6 | 7 (0,2™/ 2 )) 



(18) 



is 96% and for m = 3 it is 97%. As we will see, this 
discrepancy has a detrimental effect on the overall per- 
formance of the repeater. This could be avoided by un- 
squeezing the approximate squeezed cat states going into 
the entanglement connection. Such unsqueezing opera- 
tions may however be technically demanding, and we pre- 
fer not to include them here. We therefore consider the 
simplest situation where we directly connect the states 
generated in the first step. Alternatively, the problem 
could be mitigated by increasing m. E.g. for m > 5 we 
get an overlap of 99% with the state in Eq. (18). How- 



ever, going to such high m would also be very demanding 
in practice. 

Below we analyse the full repeater protocol, including 
entanglement swapping. Before proceeding we first ex- 
amine the performance of the connection step itself. We 
compute the output state for finite r and lossy trans- 
mission channels. Loss is modelled by fictitious beam 
splitters of transmittivity 77, such that the probability for 
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a photon to get lost on the way to the central station is 
1 — 7]. We assume that the photodetectors do not resolve 
the photon number. Details of the calculation are given 
in App. [B| First we study the output fidelity, F connect of 
the connected state with respect to |\P m ) as a function of 
the reflectivity, r of the first two beamsplitters. For this 
purpose we simulate the connection of states of the form 
in Eq. ([9| for a fixed number of iterations (m). We re- 
strict the simulations to small r since this is the relevant 
regime of the repeater. This implies that the probabil- 
ity of a successful connection is P connect ~ P c ,noioss(r)r] 
where P c ,noioss(r) is independent of the losses in the op- 
tical fibers. The results of the simulations are shown in 

Fig.m 




connect 



FIG. 4. (Colour online) The fidelity of the connected state 
with respect to the state |^ m ) plotted against the rescaled 
probability of a successful connection. We have assumed the 
input states to be of the form in Eq. (p)]). 
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FIG. 5. (Colour online) Influence of finite acceptance intervals 
in the growth on the state after connection. The fidelity of the 
connected state is w.r.t. the state |^ m ) and A is represented 
through Rgrowth- Here we have neglected losses in the optical 
fibers and R gr owth is in units of the rate at which the one- 
photon input states for the growth can be provided. 

Fig. [5] shows the same kind of behavior as Fig. [3] tak- 
ing into account that the fidelity in Fig. [3] is w.r.t. the 
squeezed one mode cat state in Eq. Q, i.e. as opposed 
to Fig. [3] the fidelity approaches unity. For optimizing 
the performance of the full repeater it is advantageous 
to have an analytical understanding of the entanglement 
generation. We have therefore fitted the graphs to func- 
tions of the form F connect = 1 — c *e d * R 9 rowth . The details 
of the fits are shown in App. [D] and Table [TJ Fig. [4] and 
Fig. [5] show that the highest rate of entanglement gen- 
eration is obtained for m = 1 but as we will see below 
we need to go to higher m for the swapping procedure to 
function. 



C. Entanglement swapping 



Fig. [4] shows that the fidelity depends linearly on 
Pconnect/v m the limit of small r. Furthermore the rate 
of the connection step for a fixed distance is more or less 
independent of m for small r. These results can be un- 
derstood by noting that the connection fails if a second 
photon is tapped off at the beam splitters. The proba- 
bility for this to happen conditioned on at least a single 
photon being tapped off is ~ Pconnect/v regardless of m. 
The second parameter to consider in the connection step 
is the vector of acceptance intervals for the growth step, 
A, which determines the fidelity of the input states with 
respect to \ipm). To determine the effect of finite accep- 
tance intervals in the growth procedure on the state after 
connection, we simulate the connection step for different 
A, taking the limit of r — » and r] — > 0. We take A to 
be the vectors giving the optimal fidelity for a given rate 
Rgrowth m Fig. [3] The result of the simulations is shown 
in Fig. [5] 



The final step of our altered repeater is to merge entan- 
gled segments via entanglement swapping. The method 
is the same as in the protocol of Ref. [27] and is illus- 
trated in Fig. [TJiii) . Two modes at the same location 
from two entangled pairs are connected on a balanced 
beam splitter and the X and P quadratures are subse- 
quently measured. Whether the swap attempt was suc- 
cessful is conditioned on the outcome of the X measure- 



ment. When swapping two states of the form (18), the 



wave functions of the states have two peaks; one at x > 
and one at x < 0. Thus, following similar arguments as 
for the growth procedure, there are two paths leading to 
outcomes in the vicinity of zero, \x\ < S. Measuring the 
plus combination there is one from the first mode having 
a positive value of x combined with a negative value from 
the second mode and vice versa. If x ~ 0, the two remain- 
ing quantum memories are projected into an entangled 
state of the form in Eq. (18) with a phase determined by 



the outcome of the P measurement. The entangled states 
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produced in the connection step are however not exactly 
of the ideal form in Eq. (18). Therefore we need to inves- 



tigate how the swapping performs with the actual states 
generated by our protocol. To this end we first identify 
the entangled state that most closely resembles the result 
of swapping after ideal growth and connection by swap- 
ping states of the form in Eq. (14). Swapping two copies 
of |^m) using the approximations (16), (17) to determine 
how the two modes gets mixed we find 



|*m)=A|0 m )|0 m )-A*|l m )|l m )+ 
C|l m )|0 m )+C*|0 m )|l m ), 



(19) 



where the coefficients depend on the measurement out- 
comes of the X and P measurements in both the current 
and previous swap levels (see App. [C] for details). This 
state contains one ebit of entanglement, and we will use 
it as our target state when evaluating the performance of 
the repeater. That is, we measure the quality of a final 
state p produced by the repeater by the fidelity 



F(p) = K* TO | / 3|* ro >r 



(20) 



The approximate form in Eq. (19) is however only 



obtained in the limit of large m. For finite m, even the 
state |^ m ), obtained in the limit of ideal growth and 
connection, will produce less than one ebit of entan- 
glement. To quantify this, we examine the dependence 
of F on the outcome of the P measurement. This 
behavior is shown in Fig. [6] where we plot F against 
the P-outcome for different values of m. For small 
values of m, there is a strong dependence. However as 
m increases the P dependence decreases because the 
states begin to resemble locally squeezed two mode 
cat states, which are insensitive to the P-outcome 
when swapped. The probability of a successful swap is 
determined by the acceptance interval 5 for the outcome 
of the X-measurement and has an upper bound of 1/2, 
which is approached for high m as in the protocol of 
Ref. [27 . Near-deterministic swapping can in principle 
be achieved following the method of Ref. [27] using 
auxiliary single-mode cat states but we will not consider 
this possibility here. 

The strong P dependence for small m in the fidelity of 
the swapped state was not seen in the original hybrid re- 
peater [27] where the the outcome of the P measurement 
merely resulted in an overall phase in the swapped state. 
As a consequence the states produced in the connection 
step of the altered repeater do not swap as well as those 
in the original repeater for the same number of iterations 
m. For long distances a large number of swap levels is 
needed. One therefore needs to go to higher m in the 
altered repeater as compared to the original repeater to 
reach a given output fidelity of the distributed state. 



0.95 




0.85 



0.95 




0.85 
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FIG. 6. (Colour online) Upper curves (left axis): The fidelity 
after entanglement swapping as a function of the P-outcome 
for (a) m — 1, (b) m — 2 and (c) m — 3, and various values 
of the X-outcome. The fidelity is F(p^) where p^ results 
from swapping two copies of |\£ m ). Lower curves (right axis): 
The corresponding probability distributions of p for each X- 
outcome. 
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IV. PERFORMANCE 

The full repeater protocol is the nested collection 
of the three steps described above i.e. growth of cat 
states, connection, and entanglement swapping. To 
quantify the performance of the repeater we use the 
fidelity F, as given in Eq. (20), and the production rate 



for the final entangled states. We set a target value of 
F > 80% and make a numerical optimization of the 
rate as a function of distance by simulating the repeater 
for different values of the control parameters at each 
step. The relevant parameters are given in Table [l[ 
We perform a full optimization over all the parameters 
in Table [l| under the constraint that the final state 
should have a minimum fidelity F > 80%. We do this 
optimization for each distance and for each value of the 
local repetition rate. For the simulation, we assume 
perfect quantum memories, perfect homodyning, and a 
SPD efficiency r] sp d = 50%. For a repeater of total length 
L and n swap levels, the distance between the stations 
is Lq = L/2 n and the transmission efficiency incurred in 
the entanglement generation step is e - L o/ 2L att ? w here 
L aU is the attenuation length of the channels. The 
total efficiency incurred is thus r] = f] sp de~ L °^ 2Latt . We 
assume L att = 20km corresponding to optical fibers 
at telecom wavelengths. The time needed for classical 
communication during entanglement generation is given 
by Lq/c, where c is the speed of light in the channels. 
We assume c = 2 • 10 5 km/s. The time required for 
local operations (measurements and memory operations) 
is assumed to be negligible compared to the classical 
communication time, such that the characteristic rates 
in the protocol are c/Lo and the source rate for the 
two-mode squeezing sources r rep . The latter is taken to 
be the repetition rate of a single two mode squeezing 
source, i.e. the rate at which down conversion is at- 
tempted in a single crystal. The optimal pair production 
probability p pa ir is found in the numerical optimization 
of the rate of the repeater for a given r rep . The effect of 
two-photon contributions in the input states is treated 
by perturbation in the pair-production probability p pa y 



as in Ref. E3 



determines the rate of the growth, 



which is the first step of the protocol and thus has a large 
effect on the overall rate of the repeater. For runtime 
reasons, we have restricted the number of growth steps 
to m < 3 and the number of swap levels to n < 4. 
For simulating both the growth and the connection step 
of the repeater we use Wigner functions to obtain the 
average output fidelity for a given set of values of the 
control parameters (see App. |A| and App. [B|. However 
this is not possible when simulating the entanglement 
swapping since the target state depends on the outcomes 
of the X and P measurements. To obtain an average 
fidelity of the entanglement step we therefore pick the 
measurement outcomes according to the probability 
distributions of X and P and calculate the fidelity of 
the resulting state. We repeat this procedure 100 times 
for each swap level and calculate the average output 



fidelity. This gives a standard deviation of the mean of 
the fidelity of about 1%. 

When performing the numerical optimization of 
the rate we calculate the fidelity of the distributed 
state and the rate on a grid of values for all the control 
parameters. The parameters affecting the performance 
of the repeater are summarized in Table [l[ In order to 
pinpoint the relevant parameter regime we use the fits 
listed in Table [I] to make an analytical approximation of 
how the fidelity depends on the different parameters. We 
use this approximation to optimize the rate using the 
method of Lagrange multipliers to find the optimal rate 
for a target fidelity of 80%. The resulting values of the 
control parameters is then used to make a grid of values 
for the numerical optimization around the analytical 
results. Finally we pick the grid point with the highest 
rate where F > 80%. The optimal rate as a function of 
distance is shown in Fig. [7] for different values of r rep . 
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FIG. 7. (Colour online) The optimal rates of the present and 
previous repeater protocols for different values of r rep . The 
protocols are optimized over the parameters listed in Table 1, 
under the constraint F > 80%. The altered repeater performs 
significantly better than the previous protocol even for r rep = 
1MHz. 

Naturally the rate of the altered repeater is very depen- 
dent on the source repetition rate. With a fast local rep- 
etition rate cat states can be grown rapidly thus remov- 
ing a time consuming step of the original repeater where 
this was done non-locally. Assuming an experimentally 
accessible repetition rate of 1MHz, the present protocol 
achieves a rate of ~ 0.08 pairs/min at L = 1000km while 
the rate of the previous protocol for the same distance 
and target fidelity is ~ 0.004 pairs/min (note error in 
[27]). The altered repeater thus gives a significant in- 



crease in the rate. For r r 



1GHz the task of storing 
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Parameter 


Description 


Effect 


Fidelity-fit 


Ppair 


Pair-production probability 
of the sources of two-mode 
squeezed vacuum states. 


Small Ppair — > low production rate 
of input states. 

Large p pa ir — > large two-photon 
component. 


F = (1 - T • Ppair)Fl + T • Ppair ' F 2 


A 


Vector of acceptance inter- 
vals in the growth proce- 
dure. 


Large acceptance intervals — > high 
growth rate. 

Small acceptance intervals — > high 
fidelity of the one-mode states state 
in |9|). 


TP — 1 r p^n.mFi growth 

1 — J- <^n,m& 


m 


Number of iterations in the 
growth step. 


High m — > low growth rate. 
Low m — )> poor swapping states. 


F = i n + j n • m 2 + k n • m 


r 


Reflectivity of the first two 
beam splitters m the connec- 
tion step. 


Large r high connection rate, 
bmall r — > high fidelity with the 
state (14). 


TP — 1 7, ( P connect \ 2 h. Pconnect 


5 


Acceptance interval in the 
swapping procedure. 


5 determines the probability of a 
successful swap and the fidelity of 
the output state. 


F e n ,me fn > mS + g n ,me hn > m \ m<2 
F = e n ,3 + fn,3 • S , m = 3 


n 


Number of swap levels. 


n determines the classical commu- 
nication time (Lo/c) between the 
stations in the elementary segments 
and hence the loss in the fibers dur- 
ing connection. 


F = 1 - lm • n 2 



TABLE I. P arameters considered in the numerical optimization of the repeater. The last column is a functional fit of how the 
fidelity (20 ) depends on the parameter when the other parameters assume their ideal values i.e. p pa ir — >• 0, A — >• 0, r — » 0, 5 — >• 0. 
r is represented through Pconnect in the fit where Pconnect is the probability of a successful connection and A is represented 
through Rgrowth- The fits for R gr0 wth, Pconnect and S are made for a specific choice of n and m. See App. [D] for details on 
the matrices containing a n , m ../i nm and the vectors containing i n --k n and l m . The expression for the fidelity's dependence on 
Ppair is calculated by perturbation in p pa ir- F\ is the fidelity with pure one-photon input states for a given set of parameters 
and F2 is the fidelity for the same set of parameters but with one of the input states being a two-photon state, r ■ 



where fi is a factor that accounts for the different acceptance probabilities for 
the repeater. 



one-photon state and a two-photon state in 



the signals in quantum memories will be challenging but 
a rate of ~ 1.5 pairs/min would in this case be reachable 
within the above assumptions. The ratio of the rate of 
the modified to that of the original repeater decreases 
as a function of the distance. This is because the states 
produced in the modified protocol are less robust to the 
swapping procedure than the states produced in the orig- 
inal repeater. When the distance increases the number 
of swap levels increase, which results in a decrease of the 
ratio of the rates for a fixed fidelity of the distributed 
state. 



V. CONCLUSION 

We have modified the quantum repeater protocol of 
Ref. [27 to improve the entanglement distribution rate. 
By interchanging the order of entanglement generation 
and growth of cat states, we have made the latter a 
local, hence faster, process, thus increasing the rate if 
local operations can be done rapidly. Furthermore, we 
have optimised the growth protocol. For entanglement 
generation, we have incorporated the method for con- 
necting cat states of Ref. [17]. We have performed a 



numerical simulation of our protocol, confirming that it 
does indeed lead to an increased rate. The final rate de- 
pends on the repetition rate of the two-mode squeezing 
sources at the base level of the protocol. For a mod- 
erate repetition rate of 1MHz, our protocol is 20 times 
faster than the repeater considered in Ref. [27], achiev- 
ing a rate of ~ 0.08 pairs/min over 1000km. This rate is 
comparable to the best proposed atomic-ensemble based 
repeaters for similar detection efficiencies (taking into ac- 
count that we have optimized for a final fidelity of 80%) 
[30] . Working with discrete variables requires SPD effi- 
ciencies of ^ 90% to obtain similar rates or complicated 
swapping procedures [15] . Much higher source repetition 
rates than 1 MHz are plausible with parametric down 
conversion in nonlinear crystals, but compatible quan- 
tum memories operating at such high frequencies may 
be very difficult to implement. For quantum repeaters of 
this kind the most feasible quantum memories are cur- 
rently those based on atomic ensembles. The high optical 
depth of a dense ensemble of cold atoms enables a strong 
coupling even for a few photons and this can provide an 
increase of the bandwidth scaling as jd : with 7 being 
the decay rate and d the optical depth. For a sufficiently 
high d the bandwidth may enable high repetition rates 
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[31]. Progress along this line was recently reported in 
Ref. [32 j, which showed memory operations with pulses 
of spectral bandwidth exceeding 1 GHz. The storage- 
and-retrieval fidelities currently achievable are far from 
the perfect case assumed in the present analysis [33 . The 
efficiency of an atomic ensemble memory can however in 
principle be made close to 100% [34]. Since the modified 
repeater does not operate with bigger cat states than the 
original repeater we do not expect different scaling of the 
two when including inefficient quantum memories. Fur- 
thermore since the modified repeater operates faster than 
the original repeater the states do not need to be stored 
for as long a time, and the effects from decoherence will 
thus be smaller [23] . Thus we expect the improvement of 
the present protocol over the previous protocol to persist 
with at least the same factor even with non-ideal mem- 
ories. It would be an interesting extension of our work 
to include non-ideal memories in the simulations, giving 
a more realistic calculation of the distribution rates but 
this is outside the scope of this article. 



Thus the state before the X measurement is 



W^ +1 (x,p,x',p') =W m [ ^+^P_+P_ , x 



W„ 



V2 ' V2 
x — x' p — p' 



V2 ' y/2 



(A4) 



Using the identity (a -\-b) 1 = Q)a s 6 2 s and collecting 
powers of x and p we can write W^+i m tne f° rm: 

wl +1 (x,p,x', p ')= J2 E 

{i,i'}=o {j,j'}=o 

EE 

k — S / — Smin 
3+j' Wax 



E E 

1=0 t'=t mu 



k - s'j \s 
3 \fj 



, (-!)• ~ s x' 



' — s' „fi+i' —k k 



l-t'J \t< 



e -^ 2 +^ 2 +P 2 +P f2 ) WijWi/j/ ^ 



(A5) 



ACKNOWLEDGMENTS 

We thank Nicolas Sangouard for useful discussions and 
we gratefully acknowledge the support of the Lundbeck 
Foundation, the Carlsberg Foundation, QUANTOP and 
the Danish National Research Foundation. 



where 



Smin = max(0, k - i), s max = min(i', fc), 
^min max(0, l-j), tmax min(j', 1). 



(A6) 



The unnormalized average output after measuring x' G 
[—A, A] is found by integrating over momentum and po- 
sition 



Appendix A: Growth of cat states 

In this appendix we describe how the growth procedure 
transforms the Wigner function of the input states. In 
general we write the Wigner function as 



2^+1 2 rn '+i 



W m (x,p) = E E W ij xl P' e ~ 



(z 2 +p 2 ) 



(Al) 



2=0 j = 



For a one photon state we have m — and the matrix 
containing Wij is 



(A2) 



The effect of the growth procedure in this representation 
is to change the size and elements of the matrix w along 
with the upper limit of the summati ons. The combina- 
tion of two states of the form in Eq. ( Al ) with variables 




x,p and x',p' on a balanced beam splitter is described 
by the transformations 



x -> -i= {x + x') , p -+ -j= (p + p') 
x' -j= [x - x') , p' - p) . 



(A3) 



OO L\ 

Jdp'Jdx'W^ +1 (x,p,x',p'). (A7) 



-oo —A 



After carrying out the integrals, we can write the unnor- 
malized state after the growth procedure as 

2 7Ti + 2 2^ + 2 

W m +i{x,p) =J2J2 «¥e~ ( * 2+p2 \ (A8) 



k=0 1=0 



with 



2Tn, + l 2 7T1 + 1 

™ ki= 2-( <+</+ ^ / )/ 2 (-i) <+</+ ^ / - 

{M'}=0 O",j'}=o 
WijWi>j>K% (A)^ J (oo), 



k-l 



(A9) 



and 



4 (*) = { 



if r < 0, 

£(^(J)/dze-"V if0<r, (A10) 

-t 



where r = i + %' — k and — ^^~*' max . We have thus 
found a simple description for the Wigner function after 
a step of the growth procedure as function of the input 
Wigner function. To fin d the state after m ste ps we start 
with the matrix in Eq. (A2) and iterate (A9) m times. 
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Appendix B: Connecting Wigner functions 

In this appendix, we describe the connection step in 
terms of Wigner functions. The state before the two 
asymmetric beam splitters with reflectivity r in the con- 
nection step is the product of the Wigner functions gen- 
erated in step one of the repeater and two vacuum states 

W rn (x,p)W rn (y,q)W vac (x f ,p f )W vac (y f ,q f ). (Bl) 



Here W m ( — , — ) has the form (Al) and W vac (x,p) = 
lg-i^+p 2 ) 

7T 

The modes described by (x,x* \p,p') are on the left 
(location A) and the modes (y^y\q^q f ) on the right (lo- 
cation B), (see Fig. ) . Before the central station it 
is only necessary to focus on the modes described by 
(x,x f ,p,p f ). Parametrising sin(# r ) = y/r, the action of 
the first beam splitter is 



x — » cos(O r )x + sm(0 r )x' , 
x' — » cos(# r )a/ — sin(0 r )x, 



(B2) 



and the corresponding transformations on the momen- 
tum variables. This results in the state 

W a i(x,x',p,p') = 

W m (cos(6 r )x + sm(6 r )x' , cos(6 r )p + sm(0 r )p') x 
W vac (cos(0 r )x' — sin(# r )x, cos(6 r )p f — sm(0 r )p). (B3) 

An additional beam splitter describing losses in the op- 
tical fibers mixes x' and p' with the vacuum mode de- 
scribed by x" and p" . We parametrise the loss by 
^/rj = cos(0j). 



x -» cos(6i)x' + sin(6>/)x // , 
x" cos(<9 z )x" - sin(^)^, 



(B4) 



and the corresponding transformations on the momen- 
tum variables. The number of photons that are lost is 
not known and consequently we trace over x" and y" . 
This produces the unnormalized state: 



W a2 (x,x' ,p,p') 



co 

dx" J dy" 



W al (x, cos(9i)x'+ sm(9i)x",p, cos(6»/)p'+ sin(0/)p") x 
W vac (coa(0i)x"- sm^iW, cos(fli)P - MW- (B5) 

The modes described by (y, y', q, q') is brought to the cen- 
tral beam splitter in the same manner producing the state 
Wb2{y, y', q, q')-The action of the central beam splitter is 



V2 ' 



y 



V2 ' 



(B6) 



and the corresponding transformations on the momen- 
tum variables. Assuming that one output mode only 



contains vacuum and the other contains anything but 
vacuum, the subsequent state is projected onto 



W vac (y',q')(l-W vac (x',p')) 



(B7) 



Consequently the state in the quantum memories after 
the connection is 



CO CO CO CO 



W ab {x,y,p,q) = ± J dx'J dy'J dp' J dq' 

— CO —CO —CO —CO 

W vac (y',q')(l-W vac (x',p'))x 

W a2 (x, (x' + y')/V2,p, (p' + q')/V2) x 

W b2 (y, (x' - y')/V2, q, (q' - p')/V2), (B8) 



where N is the normalization constant. After the inte- 
gration the resulting Wigner function can be written in 
the form 



W a b(x,y,q,p) 

^ Wstkix'iryqe 

{s,t,k,l}=0 



k I s t —x 2 —p 2 —y 2 — q 2 



(B9) 



Thi s ca n be seen by writing W a 2 and Wbi in the form 
of ( |A8| ) and evaluating the integrals using the identity 
(a + b) % — Xls=o {l) al b % ~ s as m App. A The expression 
for w st ki is rather lengthy and we shalTnot reproduce it 
here. It can, however be implemented numerically and 
thus provide an efficient description of the connection 
step. 



Appendix C: Target state of swapping 

In this appendix we outline the calculations leading to 
Eq. (19) and give the expressions for the constants that 



appear in that equation. 

We consider the swapping of two states of the type 
\^m) given in (14). The state before the swap is thus 



\^m)ab\^m)a>b> Ocfci |0 m ) |l m ) + e ± | l m ) |0 m > ) a b X 

(c 2 |0 m )| 

J-m/ ~r 6 2 -L rn }|o m )W, (ci) 

where c\ = c 2 = e\ = e 2 = 1. For generality we keep the 
coefficients named ci, c 2 and ei, e 2 since it will be impor- 
tant to consider e\^^\,2 7^ 1 in order to describe later 
swapping stages. We imagine combining modes b and a' 
on a balanced beamsplitter. Using the approximations 
(16), (17), we have the following transformations up to 
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constants of l/N^r and 1/Nf on the right-hand side 

|lm)b|lm>a' "^(| V^/im)b+|-V^/im)b)|0)a' 

+ (|v / 2/i m )a' + |-V / 2/i m )aO|0)6 (C2) 

|0 m )b|0 m ) a / 

-(1^2^)^ + 1-^2^)^)10)6 (C3) 

^|(Am+Mm)/v / 2)6|(-Mm+Mm)/v / 2)a / 

- |(-Am+Mm)/v / 2)fe|(Mm+Mm)/v / 2)a' 
+ |(Am-Mm)/\/2) b |(-/i m -/i m )/\/2) a / 

- |(-Mm-/im)/V2)6|( )/v / 2)a' 

(C4) 

~H (Am+Mm)/v / 2)6 1 (-Mm+Mm)/ V^)a' 
+ |(-Mm+Mm)/v / 2)fe|(Mm+Mm)/v / 2)a' 

- |(Am-Mm)/\/2) b |(-/i m -/i m )/\/2) a / 

- |(-Am-Mm)/v^)b|(-/im-Mm)/>/2)a'- 

(C5) 

The squeezing operators S(2)bS(2) a / should multiply the 
expressions on the right-hand side but we omit these for 
simplicity. Going to the wave function picture, assum- 
ing that cos(v / 2/imP / ) ~ cos (V^/w') , sin(v / 2/imP / ) ~ 

sin(v / 2/imP / ) 5 and tnat (! + e" 2/i -) « (1 - e~ 2A ™), we 
get the transformations 

|lm)&|lm)a' "> -2cOs( V^/wO^' 2 (C6) 

|0 m ) 6 |0 m ) a / -> 2cos(v / 2/i m p / )e" p/2 "" 2 (C7) 

|lm>6|0 m >a' -> -i(e-^^--^) +e ^K-A m )) x 

e -|(^-Am)% m (y2 /imp /) e -p /2 -^ 2 ( C8 ) 

|0 m ) 6 |l m ) ' -^i( e -^^--^) +e ^(Mm-A m )) >< 

e -|(M--Am) 2 sm ^ /imp /) e -p ,2 -^ 2 - ( C9 ) 



with 



Here is the momentum variable of mode a' and x is 
the position variable of mode b. We now perform the 
X measurement on mode b and the P measurement on 
mode a' and assume that we get outcomes p f and xq. 
The unnormalized state after the swapping is 



lm)a | Ira) 6' 

+ cic 2 C|l m ) a |0 m ) 6 ' + eie2C*|0 m ) a |l m ) 6 ' 



with 



A' = 2cos(v / 2/i m p[ ) ) 

^MraPo)- 



g 2 Am) gj^-Q 



(CIO) 

(Cll) 
(C12) 



Whenever e$ = c* this state contains one ebit of entan- 
glement and can be written as: 



A\0 m ) a \0 m } b > - A*\ 

lm)a | lm)b / H - 
C\lm)a\^m)b' + C* |0 m )a | Ira) 6' , 



A = ctc 2 2cos(V^/Wo) (C14) 
-^) 2 sin(V2/i m Po). (C15) 



Swapping two states of the form in Eq.(C13) with coef- 
ficients A\,C\ and A^^C^ respectively, it can be s hown 
within the same approximations leading to Eq. (C13) 
that the swapped state will also be of the form: 



A|0 m )|0 m/ 



A*IOIl 



m *-m 



C|l m )|0 m ) + C*|0 TO >|1 7 



(C16) 



with coefficients 



Acos(v / 2/W>o) + Bsin(V / 2/i m Po) 



C = [Ccos( V^/Wo) + Bsin(\/2/i m p / )l (C17) 



with p' again being the outcome of the P measurement. 
The coefficients A, B, C and D depend on the measure- 
ment outcome of the X measurement in the relevant swap 
and in the previous swap levels as well as t he P measure- 
ment in the previous swap levels (see Eq ( [C14] )). In the 
simulation of the repeater we replace V^MmPo 6p and 
optimize the fidelity between the swapped state and the 
target state with respect to 6 P . Note that regardless of 
this we always calculate the fidelity with a pure state 
containing one ebit of entanglement. 



Appendix D: Parameters of fidelity fits 

In this appendix we list the parameters of the fidelity 
fits shown in Table U of the article. The matrices below 
contain the constants a ny7n -h n}7n . The first entry in a 
matrix is for n = 0, m = 1 and so fourth. 






-2.19 
-9.75 






-5.39 -6.81 

-14.6 -20.1 



(C13) 



-15.6 -26.1 -39.9/ 

/0.90 0.91 0.95\ 

1.40 1.53 1.65 

2.25 3.08 3.40 

3.69 4.92 5.83 

\4.26 6.46 8.54/ 
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/0.0063 1.0 4.7 \ 

0.223 1.50 5.08 

0.460 2.59 6.26 

1.56 3.73 8.77 

V 2.02 6.68 16.1/ 

/15.0 24.2 92.0\ 

13.1 23.1 93.8 
12.3 21.6 94.5 

10.2 21.0 92.6 
\9.47 19.0 83.8/ 

/ - 
-466 -8.90-10 
0.969 -2.71-10 
1.32 -1.93-10 

\0.824 -1.74-10 

/ - - 



10" 



-7 


0.993 


-6 


0.979 


-5 


0.954 


-4 


0.905^ 



\ 



0.351 5.39 0.411 -10- 3 

0.324 5.32 0.792 -10- 3 

-0.592 4.58 3.81 -10- 3 

V-0.411 3.60 5.30 -10- 3 / 



468 0.985 

-0.969 0.951 

-0.636 0.893 

-0.260 0.799 



0.350 
0.324 
-2.45 
-4.24 



-2.80 • 10~ 3 
1.98- 10~ 3 
3.78 • 10- 3 
7.20 • 10~ 3 
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Matrices e — h show that the state's swap performance 
increases for large values of m and matrices a — d show 
that the fidelity drops as a function of r and A. The 
fact that a ny m < reflects that the fits are made for 

■^connect < ^ !• _ _ 

The numerical vectors of the constants i n -k n and l m 
are 

i = ( 1 0.938 0.811 0.618 0.413 ) 

( -0.460 -1.61 -3.14 -3.64 ) • 10" 3 
(o 0.0323 0.104 0.207 0.275 
(o.0112 0.0222 0.00542 



J 
k 
1 



Vectors i — k show that the output fidelity drops as a 
function of n and vector 1 shows that as m increases the 
states gets more robust to the swapping procedure. 
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